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Perturbative QCD calculations with heavy quarks 



1 Introduction 
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We review recent calculations of second order corrections to heavy quark decays. 
*r ■ Techniques developed in this context are applicable to many processes involving 

heavy charged particles, for example the nonabelian dipole radiation. Calculations 
involving relativistic radiating particles are also discussed. 
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Description of the higher order QCD effects in heavy quark processes is im 
' portant in view of the improving experimental precision. Ongoing efforts to 

Qh. determine the CKM matrix elements motivate a thorough theoretical descrip- 

O ' tion of B meson decays (for a review see El). These studies are made difficult 

by the complexity of the initial state, however even their conceptually simplest 
' part — perturbative 0{a 2 s ) QCD corrections to a free quark decay — were 

not known until recently. Their expected size was controversial and various 
authors used to assign widely different theoretical uncertainties in predictions 
for the semileptonic B — > X c lvi decays. 

In this talk we present a summary of our recent results on the second order 
corrections to charged fermion decays. We also discuss other applications of 
the techniques developed in those calculations. 

An exact calculation of the second order corrections to a quark (or, more 
generally, any charged particle) decay is presently not possible for technical 
reasons. On the other hand, nonejaf the standard methods of approximate 
calculations (for a short review seea) seems to be useful for those problems, 
mainly because it is not clear what could play a role of a small parameter. Here 
we will demonstrate how one can employ expansions in rather large parameters 
to obtain good accuracy of approximate results. With the example of the one- 
loop corrections to the top quark decay we will also show that in some cases 
such expansions can also give exact analytical results. 



i 



2 Computational techniques 



As an example we first consider the semileptonic decay b — > clvi. Its total 
width can be expressed as an integral over the invariant mass of the leptons 



A good estimate of 0(a 2 ) effects in T can be obtained by computing 0(a 2 ) 
corrections to dr /dq 2 at several values of q 2 . For the sake of clarity we focus 
on the point q 2 = 0; similar methods are applicable in a more general context. 

To order a 2 one has to consider three sources of corrections. Namely, there 
are virtual corrections, the one-loop corrections to a single gluon emission, and 
the radiation of two gluons. Consider the situation when the masses of the c 
quark and the b quark are very close to each other. In this case, the kinematics 
of the process simplifies considerably - the final c quark is produced almost at 
rest and moves slowly; as a result, the radiation of real gluons is suppressed. 
The calculations in the limit mt — m c <C nib can be considerably simplified. It 
is, however, more important that these calculations can be formulated in an 
algorithmic form — which shifts the burden of the calculations to the computer 
and makes it possible to obtain many terms in the expansion in 5 = 1 — m c /m{,. 

Consider for example the two-loop virtual corrections. In a general case, 
they are a function of kinematic invariants like the momentum transfer to 
the leptons, and two different masses. In the limit m c — > m,b, the virtual 
corrections simplify — they become dependent on a single parameter m,b only. 
To account for the difference nib — m c we expand the diagrams in 5. The 
coefficients of the expansion are single scale Feynman diagrams. In higher 
order terms in this expansion the propagators in those diagrams can have high 
powers. Integration by parts techniques u provide algebraic relations among 
these integrals. A solution ol those relations allows us to reduce any diagram 
to a set of "basic" integrals. □ 

Consider now the one-loop corrections to the real gluon emission. We 
essentially follow the same strategy as above. The peculiarity of this contribu- 
tion is that the Taylor expansion-.jjji S is not sufficient and one has to resort to 
the so-called eikonal expansion. Hta For the radiation of two gluons, a special 
parameterization of the phase space allows us to obtain a systematic expansion 
in S. 
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3 Decay b — > cZz/; and other applications with heavy quarks 

The techniques outlined in the preceding Section enabled us to calculate the 
0(a 2 ) corrections to the differential width dr(6 — > clvi)/dq 2 for three values 
of the lepton invariant mass QtrEn g 2 = 0, q 2 = m 2 , and g 2 = (xm, — m c ) 2 . In 
this last point the BLMEJ corrections were known previously. Using these, 
results to fit the q 2 dependence of the second order correction, we estimated H 
the second order QCD corrections to the total decay width of the b quark 
b — > clvv. 

^ = ^M\ V ^ F (4) f 1 - 1-67 ^^ - (9.8 - 1.4 ± 0.4) . 

1927T' 3 \ m b J L 71 V 7T / 

(2) 

Here F(x) = 1 — 8a; — 12x 2 In a; + 8x 3 — x and mj and m c are the pole masses 
of b and c quarks for which we use m c /m/, = 0.3. For the sake of clarity we 
separated the BLM (—9.8) and the non-BLM (1.4) parts of the second order 
corrections. We also indicated the uncertainty of our estimate of the second 
order non-BLM correction. The BLM corrections were calculated by Luke et 
al. and by Ball et al. Ii3't3 There appears to be a small discrepancy between 
these resulta^with the latter group obtaining a magnitude larger by 1.5% than 
the former. EJ This difference is smaller than the non-BLM error estimate and 
we neglect it here. 

The second order corrections in eq. (^) appear to be rather large, partic- 
ularly because of the large BLM corrections. T-his-.is related to using the pole 
quark masses in the expression for the width. Ealla It is possible to introduce 
appropriate (short-distance) masses, which reduce the magnitude of the sec- 
ond order perturbation xorrections in (^|) and improve the convergence of the 
perturbation series (seeu and references therein). 

From a more general perspective, the techniques described above provide 
an algorithm for performing the Heavy Quark Expansion of Feynman diagrams 
to order 0(a 2 ). The applicability of these methods is therefore not limited to 
the total semileptonic decay width of a heavy quark. We now discuss one 
interesting example which can be worked out using these methods. 

We consider a generalization of the QED results for the dipole radiation 
to the nonabelian case of QCD. We will later comment on why this is useful 
in practice. Let us considerl!3 a process of scattering of a color-singlet "weak" 
current J with momentum q on a heavy quark Q in the small velocity (SV) 
kinematics. For simplicity, the initial quark is assumed to rest. The initial 
Q and final Q quarks can have arbitrary masses; however, both masses must 
be large, so that the nonrelativistic expansion can be applied. The SV limit 
v = q/rriQ, \v | <C 1, is kept by adjusting q. For simplicity we consider the case 
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of equal masses, ttIq = m,Q, although nothing depends on this assumption. 
The current J must have a non- vanishing tree level nonrelativistic limit (e.g., 
Jo = QloQ or Jg = QQ), otherwise it can be arbitrary. 

Inclusive processes of the scattering of heavy quarks are described by an 
appropriate structure function W(qo, q) which is a sum of all transition prob- 
abilities induced by J into final states with momentum q and energy mq + qo- 
The optical theorem relates it to the discontinuity of the forward transition 
amplitude T(qo, q) at physical values of qo : 

T(qo, q) = d 4 xe~^ {Q\T J(x) J+(0)|Q) , W(q , q) = 2ImT(g , <f) . 

(3) 

The structure function W takes the following form in the heavy quark 
limit: 

W(u, v) = NS(u) +^^l + (w 4 ) . (4) 

At v = only the elastic peak is present. The dipole radiation is described by 
d(u>) and TV is a (velocity-dependent) factor depending on the current. 

In close analogy to QED, we define the dipole coupling ai d ^ by projecting 
Eq. (Q) on its second term: 

- Hm hm -L ^ . (5) 

iru) v^o mQ-,oo 2v z J W(uj', v)duj' 

Here to is assumed to be positive. The denominator in the last ratio eliminates 
the overall normalization of the effective nonrelativistic current. The OPE and 
factorization of the infrared effects ensures that ai d \uj) is process independent. 

Using the technique described in Section 2 we obtain the dipole coupling 
constant to second order in perturbation theory: 



3 2 A \ 6 12 



where Ca = N c for an SU(7V C ) gauge group and Po = — (N L is the 

number of flavors contributing to the running of a s ). 

An interesting feature of this coupling constant is that it determines the 
renormalization group evolution of the basic parameters jOf the heavy quark 
expansion like A(/i), /i^(/i), etc., in a scheme suggested inE3. 

In perturbation theory, the knowledge of the dipole coupling permits us to 
obtain perturbative expressions for A(/i), /^(^). These relations can be used 
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to derive a relationEZl between the pole-mass of the quark and the so-called low- 
scale kinetic mass of the heavy quarkEII to order 0(a1). The phenomenological 
importance of the low scale running quark masses is described in detail inE3. 

Similarly to the above calculation, it is possible to obtain the 0(a 2 s ) cor- 
rection to the small velocity sum rules for heavy flavor transitions—relevant for 
the determination of \V c b\ from the exclusive decays B — > D*li>iX3 



4 Decay t -> bW 

The techniques described above are useful when the quark in the final state 
can be considered non-relativistic. This obviously is not the case when the 
final quark is massless. However, such applications as the QED corrections to 
the muon decay, b — > ulvi decays, as well as the top quark decay — all belong 
to this type of decay process. In this section we discuss possible approaches 
to the calculation of the 0(a 2 B ) corrections to top quark decay t — ► bW. Wc 
will present two alternative approaches, both based on the introduction of an 
artificial small parameter. 

First, one can try to apply the same strategy as indicated above for b — ► c 
transitions to the decay t — ► bW, if one neglects the mass of the W boson. Wc 
note that this approximated, tested at the 0(a s ) corrections deviates from a 
complete result on the level of 10%. Then, one expands around configuration 
of equal t and b masses, with the expansion parameter 6 = 1 — mj/mt. In 
reality, this expansion parameter is close to one, therefore one needs to expand 
up to rather high powers of 6. To fUa 2 ) accuracy, we obtained for the top 
decay width into massless W and b: Ell 

f , 3 

T t = (i _ o.866a s - 1.69(5)a2) > (?) 

where a s is the MS QCD coupling constant evaluated at the scale mt and mt 
is the pole, mass of the top quark. More details about this calculation can be 
found incil. 

We would also like to present another approach to the same problem and 
show that in some cases introduction of an artificial expansion parameter yields 
exact analytical results after appropriate resummation is performed. Though 
we will discuss only 0(a s ) correction to top decay in what follows, a similar 
technique can be used in other ap_Dlications as well. For example, in the case 
of Z decay into massless quarks E3 this method was used to evaluate mixed 
electroweak-QCD corrections. For some of the diagrams it turned out-possible 
to sum up the expansion exactly and obtain an analytical formula IH3 which 
would have been much harder to derive directly. 
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Let us briefly describe the basic idea of the method. Consider the self- 
energy of the top quark at the value of the incoming momenta k 2 = m 2 smaller 
than the physical value of the top quark mass M. If m = M, the imaginary 
part of this self-energy gives the decay width of the top quark. For m < M 
the imaginary part is still there because of cuts through the massless b and 
W lines. The idea is to consider an expansion of the self-energy diagrams in 
S = m? /M 2 , and then take the limit S — > 1 to obtain the result for the physical 
decay width of the top quark. 





(c) (d) 

Figure 1: Diagrams whose imaginary parts determine the top width: Born level (a), and 
0{a s ) corrections (b,c,d). 



To order a s the width of the top quark is given by the diagrams in Fig. [l| 
The limit of a very heavy top is taken by neglecting b mass and replacing the 
W propagator by 

p » p v i 



P 2 

The tree level width (diagram 0(a)) is 



(8) 



r° = ^[l + £ (2-lnm 2 )]. (9) 

The calculation was performed in D = 4 — 2e dimensions and the 0(e) term is 
needed for the renormalization of the 0(a s ) correction as described below. 



6 



The one-loop QCD correction consists of the Born width multiplied by 
the wave function renormalization constant Zi plus contributions of diagrams 
(b,c,d). For Z 2 we find {C F = 4/3): 



Z 2 = C F -^M- 2e 



= C F 



16tt 2 " 
16tt 2 



M' 



-2e 



1 

+ 1 

e 



1 1 

7 ~ 2 
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V 3 V — V 



n + 1 



(10) 



If we pull out a factor mrg s e 2 Cf/(256tt SyyMyy) the contributions of the 
diagrams (b,c,d) to the width are 



r 6 = 



r, 



1 1, 2 1, „2 3 

Ye~2 lnm ~2 lnM + 4 
A/ (n + 2)(n+l)n(n-l) 



2n- 1 



E 

i=2 
gn-1 

2n 



/2(2n 2 -l) 1 



J^J V2(n + l)(n + 3) (n + 2)(n + 3) 



^ V (n 2 - l)n 2 

n=2 yy > 



n 



1 1 2 29 
1 — mm . 

4e 2 24 

After adding all four contributions and taking 6=1 we find 



(11) 



r(!) 



Cf9s 1 



13 
~6~ 



n — 1 



-3E 



2n 2 



n ^ 2 (n 2 - l)n(n + 2) 



' n(n + l) 2 (n + 2) 

G F m 3 C F a s /5 2tt 2 
V2 16tt 2 \2 3~ 



(12) 



The two sums in the above formula correspond to terms odd and even in yS, 
respectively. In obtaining them one has to shift the summation index so that 
equal powers of 6 are added together. We see that we have correctly reproduced 
the result for the one-loop correction to t — > bW. c3 



5 Summary 

We have reviewed our recent calculations of the second order corrections to 
processes involving heavy quarks. The techniques we have developed can be 
applied in many situations where the radiating particles can be considered non- 
relativistic. We have also demonstrated that some relativistic calculations can 



7 



also be performed if sufficiently many expansion terms are available. Finally, 
we have discussed some recent results on semileptonic decays of the heavy 
quark into a massless quark in the final state. 
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